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We analyze the effect of a classical noise into the entanglement dynamics between 
two particles, initially entangled, subject to continuous time quantum walks in a one- 
dimensional lattice. The noise is modeled by randomizing the transition amplitudes from 
one site to another. Both Markovian and non-Markovian environments are considered. 
For the Markov regime an exponential decay of the initial quantum correlation is found, 
while the loss of coherence of the quantum state increases monotonically with time up 
to a saturation value depending upon the degrees of freedom of the system. For the 
non-Markov regime the presence or absence of entanglement revival and entanglement 
sudden death phenomena is found or deduced depending on the peculiar characteristics 
of the noise. Our results indicate that the entanglement dynamics in the non-Markovian 
regime is affected by the persistence of the memory effects of the environment and by 
its intrinsic features. 



1. Introduction 

Quantum entanglement represents undoubtedly one of the most peculiar aspects of 
quantum mechanics as it can be viewed as the furthest departure of the quantum 
world from the classical one [T]. Furthermore, in the last years, entanglement has 
been recognized as the fundamental resource for quantum information processing 
and communication, such as quantum cryptography, teleportation and exponen- 
tial speed-up of specific computational tasks [2][3]. The single-system decoherence, 

1 



March 2, 2013 17:38 WSPC/INSTRUCTION FILE Bordone 



2 P. Bordone, F. Buscemi, C. Benedetti 

due to the environment, can make the entanglement disappear on the short-range 
time scale, making the quantum parallelism, essential for quantum computation, 
ineffective [4]. On the other hand, environment can even resume entanglement or 
preserve it. Thus, it is still very important to analyze the effect of the various kinds 
of environmental noise on the entanglement dynamics in realistic quantum systems 
which can exhibit peculiar phenomena, such as entanglement sudden death (ESD), 
and entanglement revival (ER) [5 . 

ESD, namely the disappearance of the quantum correlations between two quan- 
tum systems at finite time in spite of an exponential decay of the system coher- 
ence, appears to be peculiar of the Markov environments, namely environments 
with short, or rather instantaneous, self-correlations [BJ. Instead, ER, that is the 
rise up of quantum correlations in a bipartite system after a finite time period 
when they completely disappear, occurs in presence of non-Markovian noise [7], 
which represents the dominant source of decoherence in solid-state systems. The 
presence of environments with memory [7] is crucial for the non-monotonic time de- 
pendence of the entanglement which cannot be ascribed to interaction between the 
two sub-systems. Indeed, ER is found both for couples of quantum systems inter- 
acting directly or indirectly in a common quantum reservoir and for noninteracting 
quantum systems in independent non-Markovian quantum environments. 

The effect of the quantum noise on the entanglement dynamics between two 
quantum systems has been interpreted in terms of the transfer of the correlations 
back and forth from the systems themselves to the environment. This is due to the 
back-action of the systems on the environment. However, recent works showed that 
the non-monotonic time dependence of the amount of quantum correlations may 
occur in two-qubit systems under the local action of a system-unaffected environ- 
ments, such as classical random external potentials [HE] - Due to the classical nature 
of the noise, in this case no back-action-induced correlations can be transferred from 
the system to the environment. Thus, the occurrence of ER in a non-Markovian 
classical environment is in contrast to the well-established interpretation of revivals 
in terms of system-environment quantum back action and rises the fundamental 
question of how one could explain the effect of classical noise on entanglement 
dynamics in bipartite systems. 

In this work, we intend to analyze the role played by a classical noisy environ- 
ment into the entanglement dynamics between two quantum particles in a simple 
model, which not only allows us to relate the intrinsic features of the two-particle 
dynamics to the revival of the correlations, but can also be of great relevance in 
various physical phenomena. The diffusion of two quantum particles along a one- 
dimensional lattice in the presence of static and dynamic disorder, mimicked by 
a noisy potential, represents a good candidate for our theoretical investigation. In 
fact, the spreading of a quantum-mechanical wave packet in a tight-binding model 
with a stochastic potential is able to describe electronic transport in solid-state 
systems, particle diffusion in molecular crystals, diffusion of excitations, photon 
propagation in coupled waveguide lattices. Furthermore, it gives the possibility to 
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connect the transition from quantum to classical behavior of the two particles to the 
features of the environmental noise (correlation time), as suggested by investiga- 
tions on the diffusion of single-particle wavepackets [10] . This connection could be 
viewed as a valid guideline to study in details, and, possibly, to clarify how Markov 
and non-Markov classical noise affect the occurrence of ESD and ER. 

In our approach, the dynamics of the two-particles subject to a random po- 
tential will be faced numerically. The effect of noise will be simulated both by 
assuming random and time-independent amplitudes for the transitions from one 
site to another and by using a sort of dynamic disorder, expressed by a potential 
term describing a random telegraph noise. 

The paper is organized as follows. The physical model adopted is described in 
Sec. 2. In Sec. 3 are presented the results for the various cases considered: static and 
dynamic noise, Markovian and non-Markovian environments. Conclusions, com- 
ments and future perspectives are given in Sec. 4. 

2. Physical Model 

A model of two distinguishable particles continuous-time quantum walks (CTQW) 
in ID graphs is solved by means of numerical techniques and the time evolution 
of the wavefunction describing the system is then used to estimate the degree of 
quantum correlations. 

We study CTQW on a ID ring lattices of TV sites (with N even) with periodic 
boundary conditions. In agreement with previous single-particle investigations 
the topology of the graphs here considered is simple, in the sense that each node is 
connected to its two first neighbors, still it is good enough to describe experimental 
implementations of CTQW, such as two-photon transport in an array of waveguide 
lattices where non-classical correlations appear. Therefore our model represents a 
valid tool to analyze the amount of quantum correlations appearing in two-particle 
CTQW. 

The two-particle Hamiltonian describing the dynamical evolution of the system 
is given by 

H = H A + H B , (1) 
where 'H A ^ * s * ne single-particle Hamiltonian acting on the particle A(B): 

N 

WV)=EK (B) ti)MB)(j\ - cf B) (t)(\j + l) A(B) (j\ + \j) A(B) (j + l\)}. (2) 
J=l 

The kets \j)a(b) indicate the quantum states describing the particle A(B) localized 
in the j node and form a complete, orthonormalized basis set, which span the whole 
accessible Hilbert space. (3f" B ^ is the on-site energy, and cf" B ^ (t) is the tunneling 
amplitude between nearest neighbors 12 . Due to the periodic boundary conditions, 
here we assume that the site N + 1 coincides with site 1. 
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To evaluate the dynamics of the entanglement between the two particles and 
of the decoherence induced by the noisy environment on the two-particle state, we 
need to calculate the density matrix of the system as a function of time. To this 
aim we need to estimate the evolution operator of the system which is a function 
of the eigenvectors and eigenvalues of the Hamiltonian which, in turns, depends 
on the parameters j3^ B ^ and it). Therefore, once a suitable choice of such 

parameters is done, we diagonalize numerically the Hamiltonian for specific times 
and get the corresponding density matrix values. 

In order to describe the effect of the environment on the two-particle quantum 
state, one can introduce noise either through the coupling constant j3^ B ^ or through 

c^ B \t). In agreement with previous investigations [321[T3], we adopt the latter. 
Once noise has been included in the model, the state of the system is represented 
by a density matrix {pit)) which is the result of an average over a number of density 
matrices each one obtained with a specific choice of the cf^ B \tYs. In particular, to 
estimate (p{t)), first we evaluate a single numerical run by generating a sequence 
of transitions with random amplitudes and then by solving exactly the two-particle 
Schrodinger equation for that given sequence. The final numerical simulation is 
found by producing a large number of runs, each with a different sequence, and 
then calculating the average density matrix over all the runs. 

Decoherence, which is an evaluator of the degree of entanglement between the 
two-particle system and the noisy environment, is here estimated by means of the 
von Neumann entropy |14j of (p{t)): 

e vN (t) = -Tr{(p(t))\n(p(t))} . (3) 

Note that such calculation implies the numerical diagonalization of the N x N 
matrix (p), that becomes more and more demanding with increasing the number 
of sites of the graph. 

The particle-particle entanglement is quantified in terms of the negativity |15j : 

AT=^|A i ((p^)))|-l, (4) 

i 

with Aj ((p TA(B) )) the eigenvalues of (p TA ^ B >) which is the partial transpose related 
to the subsystems A{B) of the total density matrix (p). Negativity has been shown 
to be a valid measure of the degree of quantum correlations only between two 
quantum sub-systems each having two degrees of freedom, namely qubits. In fact 
no measure discriminating separable from entangled states for a mixed state with 
more than two degrees of freedom is known. As a matter of fact, non-zero negativ- 
ity, for the general case of N degrees of freedom, is just a sufficient condition for 
entanglement [IB] . 

3. Numerical results 

In this section we present the numerical results obtained for the case of N = 4, 6, 8 
sites for two different ways of modeling classical noise, namely through static and 
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dynamic disorder. In particular the former is simulated by assuming random and 
time-independent values for each , while the latter is modeled by selecting the 
(i) according to a random telegraph signal. 



3.1. Static noise 

In the static disorder case, the parameters c^^'s are assumed to be random vari- 
ables following the flat probability distribution P(c) = ^ for |c — c | < y an d 
otherwise [13j . Co is the mean value of the distribution and A is a measure of the 
disorder of the environment, in fact, when A goes to zero, the noise effect van- 
ishes. The autocorrelation function of the cf s, given by (5c(t)Sc(0)) = ^j, is 
time independent and this implies that its power spectrum is proportional to a 
(^-function centered on zero frequency. As a consequence, this kind of noise has a 
characteristic time which is always much longer than the characteristic time of the 
system-environment coupling. Therefore, the static disorder can be considered as 
representative of a non-Markovian noise. 




— N=4 
-- N=6 



t (arb.unit) t (arb.unit) 



Fig. 1. Left panel: Entanglement (solid line) and decoherence (dashed line) as a function of time 
for N = 2. Right panel: decoherence time evolution for N = 4 (solid line) and N = 6 (dashed 
line). 



In FigQ] the time evolution of decoherence and entanglement are shown at 
N = 2,4,6, for the input maximally entangled state \\p(t — 0)) = Ka*b)- 
The entanglement, here expressed in terms of negativity, as above explained is 
calculated for the N — 2 case only. From the left panel of FigO] we see that the 
entanglement exhibits ESD and ER phenomena [8], corresponding to an oscillat- 
ing behavior of the decoherence curve (here oscillations are much less pronounced 
than for the entanglement case, because of the logarithmic dependence displayed 
in Eq.Q). By increasing the number of sites (right panel), the saturation value 
of the decoherence increases because, as a consequence of the environment action, 
the system is described by a statistical mixture involving a larger number of terms, 
and, for this kind of noise, we observe that these oscillations tend to disappear. 
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3.2. Dynamic noise 

In the dynamic disorder case the parameters c^ B ' , s are set equal to vQj (t) 
where Q(t) is described by a stochastic process of the telegraph noise type, 
which flips randomly between the two values 1 and -1 at a rate 7, and v is the 
particle-environment coupling constant. The autocorrelation function of Q{t), as 
well known, is given by (5Q(t)5Q(0)) — e~ 2lt with the Lorentzian power spectrum 
(SQSQ)^ — 47/ (lj 2 + 47 2 ). The relative length of the two time scales r e ~ I/7 
(characteristic time of the environment) and t c ~ 1/v (characteristic time of the 
particle-environment coupling) are critical for decoherence and entanglement. In- 
deed, r e < t c is identified as the Markovian regime, while r e > t c is identified as 
the non-Markovian regime [8] . 

3.2.1. Markov case 

In FigJ5]the results for the case r e /r c = 1/5 are presented, for the same initial state 
|V"(0)) used in the previous section. Unlike the case of static disorder, because of 
the short time correlation of the environmental noise, the degree of entanglement 
between the two-particles decays continuously and no ESD and ER phenomena 
are observed. Decoherence is a monotonic function of time, and reaches asymptotic 
values larger and larger as the number of nodes increase. 




1 2 3 4 5 6 2 4 6 8 10 12 14 
t (arb. unit) t (arb. unit) 



Fig. 2. Left panel: Entanglement (solid line) and decoherence (dashed line) as a function of time 
for N = 2. Right panel: decoherence time evolution for N = 4 (solid line) and N = 6 (dashed 
line). Here the Markov regime is obtained by setting t s /t c = 1/5 



3.2.2. Non-Markov case 

The non-Markovian behavior has been simulated by setting the ratio r e /r c = 10. 
The results are displayed in Fig|3] For N = 2 (left panel), the entanglement is a 
damped oscillating function of time, thus showing ESD and ER phenomena. Since 
the envelope is exponential we can assume that revival happens an infinite number 
of times, in agreement with previous findings [8]. Decoherence, unlike the Markovian 
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case, presents evident oscillations before reaching the saturation value. In particular 
we note that the local minima of decoherence correspond to the local maxima of 
entanglement, thus suggesting a strict connection between the two quantities. In 
fact, in our model, we can attribute the disentanglement of the two particles to 
the movement of the initial state |^(0)) towards mixed state. Such a relation can 
reasonably be thought to hold also for graphs with N > 2, where no rigorous eval- 
uations of the entanglement can be performed. Therefore, the oscillating behavior 
of decoherence observed for N = 4, 6 (right panel of FigEJ) can be considered as an 
indicator of the presence of ER. 




Fig. 3. Left panel: Entanglement (solid line) and decoherence (dashed line) as a function of time 
for N = 2. Right panel: decoherence time evolution for N = 4 (solid line) and N = 6 (dashed 
line). Here the non-Markov regime is obtained by setting r e /r c = 10. 



4. Conclusions 

We have investigated the effect of a noisy environment on the time evolution of the 
entanglement between two-quantum particles and of the loss of coherence of the two- 
particle system for CTQW in a one-dimensional ring lattice with periodic boundary 
conditions. The two-particle initial state has always been chosen as a maximally 
entangled one. The environmental noise has been simulated introducing both static 
and dynamic disorder in the lattice. To this aim the transition amplitudes from one 
site to another have been assumed random and time independent (static case) or 
subject to a random telegraph signal (dynamic case). 

In the dynamic case, the persistence of the memory effect of the environment is 
related to the switching rate of the random telegraph signal. Therefore two distinct 
regimes can be identified according to the ratio between the correlation time of 
the environment and the one characteristic of the system-environment interaction. 
If such a ratio is smaller than 1, Markovian behavior can be assumed, on the 
contrary, values larger then 1 identify non-Markovian regime. The static disorder 
represents an environment whose memory effect can not be neglected at any finite 
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time, therefore the characteristic time of its correlations is always larger than the 
one of the environment-system coupling. Thus we classify this kind of noise as non- 
Markovian. It should be noticed that such model is different from the one adopted 
in the literature |17H19| where the non-Markovian noise is commonly represented 
only in terms of random telegraph signal. 

For the Markovian case we found an exponential decay of the entanglement, 
while decoherence increases monotonically with time up to a saturation value, de- 
pending upon the number of lattice sites. No peculiar phenomena, like ESD or ER, 
are observed, in agreement with previous findings |S]. 

For the non-Markovian behavior we find that the dynamics of entanglement and 
decoherence exhibits non-monotonic time dependence. In particular ESD and ER 
are detected [THS] for N = 2, and can even be supposed to be present for N > 2 
under some specific conditions (dynamic noise). Indeed our approach suggests that 
different kinds of noise leads to somehow different results. For static disorder both 
entanglement and decoherence show less pronounced oscillations with respect to 
the non-Markovian random telegraph noise for N — 2, while for N > 2 decoherence 
increase practically in monotonic way, thus suggesting that ESD and ER phenomena 
are not present, unlike what is found for the dynamic noise case. 

In conclusion, our analysis points out that the entanglement of a quantum sys- 
tem coupled to a classical source of noise is affected by the persistence of the mem- 
ory effect of the environment and by the intrinsic features of noise itself. Therefore 
we believe that the extension of the present investigation to the ubiquitous case 
of l/f noise could be of great interest to get a better insight into the effect of a 
classical noise on a quantum system [20] . Furthermore, in order to get a deeper 
understanding of how a classical environment affects the appearance of quantum 
correlations, the study of more complex systems, including many-body effects and 
a higher number of degrees of freedom, is required. In this view, the extension of 
the approach here adopted to model many particles interacting with each other is 
not straightforward and will need the development of more sofisticated and efficient 
numerical tools. 

Acknowledgements 

The authors would like to thank M.G.A. Paris for fruitful discussions. 



References 

[1] E. Schrodinger, Naturwissenschaften 23 (1935) 807. 

[2] F. Buscemi, P. Bordone, and A. Bertoni, Phys. Rev. B 81, (2010) 045312. 

[3] F. Buscemi, Phys. Rev. A 83, (2011) 012302. 

[4] F. Buscemi, P. Bordone, and A. Bertoni, Phy. Rev. A 73, (2006) 052312. 

[5] T. Yu and J.H. Eberly, Science 323, (2009) 598. 

[6] T. Yu and J.H. Eberly, Phy. Rev. Lett. 93, (2004) 140404. 

[7] B. Bellomo, R. Lo Franco and G. Compagno, Phy. Rev. Lett. 99, 160502 (2007). 

[8] D. Zhou, A. Lang, and R. Joynt, Quant. Inf. Comp. 9, (2010) 727. 



March 2, 2013 17:38 WSPC/INSTRUCTION FILE Bordone 



Effect of Markov and non-Markov Classical Noise on Entanglement Dynamics 9 

[9] R. Lo Franco, E. Andersson, B. Bellomo and G. Compagno, Phys. Rev. A 85, (2012) 
032318. 

[10] A. Amir, Y. Lahini and H.B. Perets, Phy. Rev. E 79, (2009) 050105(R). 
[11] O. Miilken and A. Blumen, Phys. Rev. E 71, (2005) 036128. 

[12] Y. Lahini, Y. Bromberg, D.N. Christodoulides and Y. Silberberg, Phys. Rev. Lett. 
105, (2010) 163905. 

[13] C. Thompson, G. Vemuri, and G.S. Agarwal, Phys. Rev. A 82, (2010) 053805. 

[14] A. Peres, Quantum Theory: Concepts and Methods (Kluwer Academic Publishers, 

Dordrecht, 1998). 
[15] G. Vidal, R.F. Werner, Phys. Rev. A 65, (2002) 032314. 
[16] C. Sabm and G. Garcia- Alcaine, Eur. Phys. J. D 48, (2009) 435. 
[17] L. Mazzola, J. Piilo, and S. Maniscalco, Phys. Rev. Lett. 104, (2010) 200401. 
[18] L. Mazzola, S. Maniscalco, J. Piilo, K.-A. Suominen, and B. M. Garraway, Phys. Rev. 

A 79, (2009) 042302. 

[19] F. F. Fanchini, T. Werlang, C. A. Brasil, L. G. E. Arruda, and A. O. Caldeira, Phys. 

Rev. A 81, (2010) 052107. 
[20] A. De, A. Lang, D. Zhou, and R. Joynt, Phys. Rev. A 83, (2011) 042331. 



